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■ Abstract. In this paper, we investigate some consequences of multiple-field inflation 

, for the cosmic microwave background radiation (CMB). We derive expressions for 

' the amplitudes, the spectral indices and the derivatives of the indices of the CMB 

power spectrum in the context of a very general multiple-field theory of slow-roll 
^vq \ inflation, where the field metric can be non-trivial. Both scalar (adiabatic, isocurvature 

J> ■ and mixing) and tensor perturbations are treated and the differences with single-field 

inflation are discussed. From these expressions, several relations are derived that can 
be used to determine the importance of multiple-field effects observationally from the 
[~^. ■ CMB. We also study the evolution of the total entropy perturbation during radiation 

and matter domination and the influence of this on the isocurvature spectral quantities. 

m 
o 

_C2 - 

PACS numbers: 98.80.Cq, 98.70.Vc 

6 
& 

in 
c3 



1. Introduction 

One of the greatest successes of the concept of inflation [TJ El HI IS] is that it can give an 
explanation for the existence of small density fluctuations in an otherwise homogeneous 
universe. These small fluctuations, which are the gravitational seeds for the formation 
of large-scale structures, are observed in the cosmic microwave background radiation 
(CMB). Of course this works two ways: the observed amplitudes and slopes of the 
fluctuation spectra give us some observational constraints on the otherwise rather elusive 
parameters in inflation models, and thus on the parameters in the underlying high-energy 
theories. Hence it is important to have expressions for these spectral quantities in terms 
of the parameters of the inflation models. The most accurate observations of the CMB 
to date are those from the Wilkinson Microwave Anisotropy Probe (WMAP) mission 
[H| (see 0012!). 

This paper is a sequel to our paper jTUj . In that paper a general theory of inflation, 
with an arbitrary number of scalar fields that take values on a possibly non-trivial field 
manifold, was considered (motivated by (string-inspired) supergravity theories where the 
Kaliler potential leads to non- minimal kinetic terms). We treated both the background, 
paying special attention to the concept of slow roll, and the scalar perturbations. After 
deriving expressions valid to first order in slow roll for the perturbation quantities at 
the end of inflation, we also considered what happens after inflation. This finally led to 
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expressions for the correlators of the adiabatic and isocurvature scalar perturbations, 
as well as of the mixing between them, at the time of recombination (basically these 
are the scalar amplitudes of the CMB power spectrum). An example of a quadratic 
potential was included to illustrate the theory. 

In this paper, these results are extended in basically two ways. In the first place, 
expressions for the spectral indices n and their derivatives dn/d Ink are given next to 
those for the amplitudes, and results for the tensor perturbations are added as well. 
Using these results we derive an expression that can be used to observationally test the 
importance of multiple-field effects in general, as well as various consistency relations for 
the one- and two-field cases. Secondly, the behaviour of the total entropy perturbation 
after inflation is studied to investigate under which conditions the usual assumption of 
a constant isocurvature perturbation is justified and what happens to the isocurvature 
spectral quantities if this is not the case. 

A large number of references to other papers on inflationary perturbation theory was 
given in the introduction of ^Oj- More recently, a second-order treatment of the single- 
field case was given in [II] , while the authors of G3] used non-slow-roll techniques 
in the single-field case. Many aspects of the two-field case were recently discussed in 
|14[ ITS'] \W\ ITTj (especially the last of these contains a very large number of references 
to earlier papers). An extensive quantum treatment, including backreaction, of O(N)- 
symmetric multiple-field inflation models was given in ^H] • Previous work on consistency 
relations for two-field inflation can be found in [THJ I2D1 EU ED EH] 

The outline of the paper is as follows. Section 121 deals with the general theory of 
multiple-field inflation. In three subsections the background, scalar perturbations and 
tensor perturbations are discussed. Most of this (except for the tensor perturbations) 
was derived in [10] . so that this section is kept brief, only setting up the notation and 
giving the results that are used in the rest of the paper. Derivations are only given for 
those results that were not present in [10] . In section E2 (plus the appendix), expressions 
for the scalar (adiabatic, isocurvature and mixing) and tensor amplitudes, spectral 
indices, and the derivatives of the spectral indices are derived. The results are discussed 
with special attention to the differences with single-field inflation. In subsection 13.31 a 
number of important relations between these quantities are derived that are relevant for 
observationally testing the significance of multiple-field effects. Section H] is devoted to 
the total entropy perturbation, which enters as a source term in the equation of motion 
for the gravitational potential after inflation in a system with multiple components 
originating from multiple fields during inflation, and thus determines the isocurvature 
perturbation. An equation for its time derivative is derived and the resulting time 
dependence is discussed. Finally, section [5] contains the conclusions. Some of the results 
in this paper have appeared in [22] . 
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2. Multiple-field inflation 

After defining the general setup of the theory, the results for all the relevant perturbation 
quantities at the time of recombination are given in this section. For the derivation of 
the scalar results as well as for more details about the background quantities the reader 
is referred to our previous paper [10j the derivation of the tensor results is given here. 
Readers interested in seeing more intermediate steps in the calculations might also want 
to take a look at [22]- 

2.1. Background 

We consider the situation where the matter content of the universe consists of an 
arbitrary number of real scalar fields a , which are represented as a vector <p of scalar 
field components. These fields are the coordinates on a possibly non-trivial field manifold 
with field metric G. The field is separated into a classical homogeneous background part 
and a small quantum perturbation part: full (?7, x) = 4>{rj) +54>(r), x). All equations are 
linearized with respect to the perturbations. The associated Lagrangean with potential 
V is 

C = V=9 ("^ • W - V(4>)) = [~ l -g^d^ T GdA - V{ct>)^ , (1) 

where T denotes the transpose and g^ u is the spacetime metric. Because of the non- 
trivial field metric we have to define covariant derivatives instead of normal ones: 
V b A a = A a b + Tl c A c and V^A a = d^A a + T^d^A" for the field derivative and the 
spacetime derivative of a vector A a in field space. 

For the metric part of the universe we make the following definition: 

^-^M" AM* 'MsO- 

where we have applied the longitudinal and vector gauges and used the (^-component 
{i 7^ j) of the Einstein equation to set \I/ = <£> in the scalar part (see e.g. [221 121] )• The 
first term represents the homogeneous background (a flat Robertson- Walker metric with 
scale factor a(r])), the second term the scalar perturbation (with $(77, x) the gravitational 
potential), the third term the vector perturbation (with S(r], x) a divergenceless vector), 
and the fourth term the tensor perturbation (with hij(r],x) a symmetric transverse 
traceless tensor, i.e. hij = hji, h 1 ^ = and h\ = 0). As proved in [213 [21] scalar, vector 
and tensor perturbations decouple to linear order and can be considered separately. We 
make use of three different time coordinates: comoving time t, conformal time 77 and 
the number of e-folds N, defined as 

dt = adr], dN = Hdt = Hd-n, H = -, H = - = aH, (3) 

a a 

where H and 7i are the Hubble parameters associated with t and 77, respectively, and a 
dot (prime) denotes a derivative with respect to comoving (conformal) time. 
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For giving a physical interpretation of the various scalar field components, it is 
very useful to define a basis on the field manifold. (It is also an important ingredient 
of the calculations, especially quantization, as explained in [10J, where we introduced 
this basis.) We define the first basis vector e.\ as the direction of the field velocity: 
e.\ = 0/|0|. Next, e.2 is the direction of that part of the field acceleration T> t <p that 
is perpendicular to e 1; and this is extended to higher-order derivatives to define the 
other basis vectors. Once this basis is defined we can use projection operators to split 
vectors into their physical components. The most important ones are P" = e 1 ef G and 
P- 1 = 1 — P" that make it possible to distinguish between effectively single-field (ei) 
and truly multiple-field effects. The basis vectors are not constant in time; their time 
derivatives are given by 

/An+l) ~(n) \ 

Ti tt I ' n+1 ' n \ / a \ 

V ^ = H ~(nj~ e n+l ~ 7(^=1) en-! , (4) 

V Vn Vn-1 J 

where fj n = e n ■ r)^ n \ with f)^ defined below. 

The multiple-field generalization of the single-field slow-roll approximation was 
treated in ^Hj as well. We defined the following slow-roll functions: 

and the short-hand notation f) = f)^ and £ = f]^. Here k is the inverse reduced 
Planck mass: k, 2 = 8ttG = 8n/Mp. Using the basis discussed above we can take the 
components rjr = e.\ ■ fj and fj 1 - = e 2 fi (by construction there are no other components), 
and similarly for although there one has not just a single perpendicular direction, but 
two (£2 and £3). Using these definitions the equations of motion for the background 
quantities are: 

K_ ' 1 ^ 1/2 

y/3 v ~ , - .... 3, 

where we have switched to using the number of e-folds iV as time coordinate because it 
simplifies the equation. These equations are still exact; the slow-roll functions have only 
been used as short-hand notation. However, the next step is the slow-roll approximation: 
we assume that the slow-roll functions are small so that we can set up a series expansion 
in terms of them, (e and 77 are assumed to be first-order quantities, while f]^ is of order 
n — 1.) Note that the time derivatives of the slow-roll functions are of one order higher 
in slow roll: 

e jN = 2e(e + 77"), V N f)^ = f) in+ ^ + ((n - l)e - fj ll W n) , 

(fjW) iN = f II + (r^) 2 + (e - r)")/]", (^) >JV = 6 + (e - 2^)rj\ (7) 

2.2. Scalar perturbations 

The derivation of expressions for the scalar perturbations was the main topic of 
[Hi] . Because the equations of motion for the gravitational potential and the field 
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perturbations do not decouple in the multiple-field case, we cannot ignore the latter 
as one can in the single-field case. The scalar-field perturbations are redefined and 
quantized as 

ty + i^'V q = Qa^ + Q*a, (8) 

with a) and a constant creation and annihilation operator vectors and Q{rf) a matrix 
function that satisfies a classical equation of motion. Here we switched to working in 
the explicit basis: non-bold versions of a vector or matrix represent those quantities in 
this basis, for example, q T = (qi, q2, ■ ■ ■) = (ei • q, e 2 • q, . . .). This means that the non- 
bold basis vectors are simple constant vectors: ej = (1, 0, 0, 0, . . .), e 2 = (0, 1, 0, 0, . . .), 
etc. In the rest of this section we will denote expressions like the one for q in (JHJ) as 
q = Qa) + ex., with the c.c. meaning complex conjugate. 

The results for the adiabatic and isocurvature parts (see below for a definition as 
well as for some remarks about assumptions made in the derivation) of the gravitational 
potential in terms of Fourier modes at the time of recombination, valid to first order in 
slow roll for modes with k -C 7^ rec (super-horizon modes), are 

3 K> II 

$ fc ad(trec) = 5^^= K + U Pe) ^4 + «, (9) 



$ fcis o(t rec ) = ~ l^T2^ V ^ E ^l + CC (10) 

Here we have defined E H = (1 — e n )l + (2 — 7 — In 2)5-^, with 7 the Euler constant, and:): 

M 2 
3H 2 



M 2 

5 = ?1 - ^EH + 2 ^ e i e L M 2 = G -1 V T W - R(0, 0), (11) 



l$ e = 2y/hl [ E dt>H^^e T 2 QQ n \ V e T ee e T 2 Q e Q H \ (12) 

The subscript Ti means that the quantity has to be evaluated at the time t H during 
inflation when 7i = k (often called the time of horizon crossing). The quantity Qn is 
given by Q-h = E-ft/\/2k. Note that the slow-roll approximation was only used during 
a small transition region around the time t^. The time t e is the end of inflation. The 
vectors Up e and V e have no components in the e\ direction. The matrix R is the 
curvature tensor on the field manifold and [R(0, 0)] a ^ = R^cdr '0 C '■ 

The terms adiabatic and isocurvature relate to the evolution after inflation. During 
radiation and matter domination the equation of motion for super-horizon modes of $, 

$ + AH(1 + f c 2 )$ + k 2 ( P c 2 - p)$ = -2k 2 ( P c 2 - p)S, (13) 

has an inhomogeneous source term proportional to the total entropy perturbation S: 

s l*Z^£ (14) 

I For normalization reasons, we have removed a factor i from the definition of V e as compared with 
our original definition in |10j. 
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with p and p the total pressure and energy density and the sound velocity c 2 s = p/p. The 
adiabatic perturbation is the homogeneous solution for $ of this equation of motion, 
while the isocurvature perturbation is the particular solution generated by the source 
term, with the initial conditions that it is zero and has zero derivative at the beginning of 
the radiation-dominated era. If S is constant on super-horizon scales, we have the simple 
solution <£> iso — ~S during the matter era, which was used in the derivation of equation 
(jlOj). To perform the matching of these solutions with the ones at the end of inflation 
(to determine the constants of integration) we assumed an immediate transition to a 
radiation-dominated universe at the end of inflation, ignoring (p)reheating. Especially 
for the isocurvature perturbations this is expected to be a crude approximation (since 
V e depends crucially on what happens at the end of inflation, see ([12))). which needs 
further improvement, but the treatment of the perturbations during a more realistic 
transition at the end of inflation as well as during an epoch of (p)reheating is still under 
investigation^ 

The behaviour of 5* during radiation and matter domination is investigated in 
section 01 to see if the assumption of constancy on super-horizon scales is justified. 
However, let us discuss here what happens with (fTT)|) in the case that S is not constant, 
but given by S(t) = f(t)S e . Here f(t) encodes both the possible time dependence of 5* 
during radiation and matter domination and a possible correction factor to the simple 
matching condition described above, and S e is the total entropy perturbation at the end 
of inflation, given by [TUl 122] 

(15) 

A particular solution for the inhomogeneous equation (|T3|) is now given by $ par t(^) = 
—2F(t)S e , where (w = p/p) 

m . 1 f itll^L A. (•) f df m JL (L±) . (16) 

w a J I + \w df \HJ J JK ' dt" V 1 + w J K ' 

Adding the homogeneous solution $hom = C{H/a) + D{H/a) f'dt'(l + w)a, with 
appropriate constants C and D to satisfy the initial conditions, we finally obtain the new 
result for $i SO (^rec)- It turns out that we can represent it by exactly the same expression 
as in equation (JTHJ) , if we change the definition of V e to 

V T = 

e 

with t* the beginning of the radiation-dominated era. Note that V e has now lost 
the physical interpretation of a quantity completely determined by inflation, but 
mathematically it is very convenient to work with an expression for $i SO (^rec) i n which 
only V e changes, as opposed to changing the form of (fTT)|) while keeping the definition of 

§ Some study on the effects of preheating on the scalar perturbations has been done for specific models, 
see 12ZI EHl 1210 EHH RH and references therein, but different authors do not yet agree regarding 
the conclusions. 
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e n^. — z\\ — e 2<2 e Qw> ( 17 ) 



e e + V< 
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V e fixed. In the case that f(t) = 1 for all t, one finds F(t) = 1 for all t, and (|17|) reduces 
to (H2J), as it should. 

2.3. Tensor perturbations 

Having reviewed the scalar perturbation results from ^H] i n the previous subsection, 
we now turn to tensor perturbations. Scalar-field perturbations cannot, by definition, 
generate tensor perturbations in the metric. However, the two tensor degrees of freedom 
of the metric are the only physical ones (representing the two polarizations of the 
graviton) and they do not need to be generated by a matter source. || Because the 
scalar fields do not generate the tensor perturbations, there is no difference between the 
treatment of these perturbations in multiple-field or in single-field inflation. Inflation 
enters only by way of the background quantities. Hence the results derived in this 
subsection are not new (see e.g. [3311231 EI] and references therein), but we derive them 
here using the methods and definitions of 

The (ij) -component of the Einstein equation gives the well-known equation of 
motion for the tensor perturbation: h'L + — A/iy = 0. Because h^ is symmetric, 

transverse and traceless, it has only two independent components, which are represented 
by two constant polarization tensors e^, A = 1,2, normalized as ef^ B = 5 AB and 
satisfying the same three properties as hij. The tensor perturbation is written as 
hij(r], x) = Y^A=i(^ K '/ a )' l l ! A{v^ x ) e tji where the factor 2k/ a has been taken out to simplify 
the equation of motion and to obtain the correct normalization of the Lagrangean. After 
switching to Fourier modes we find the usual equation of motion for the mode functions 



This equation is similar to equation (34) in JU] for the redefined gravitational potential 
Uk, but without an inhomogeneous multiple-field term. An important difference between 
ipA and u is that the ipA represent two physical degrees of freedom and can be quantized 
directly, while u is not a physical degree of freedom, and had to be quantized indirectly 
by means of the scalar field degrees of freedom q. There is no coupling between the two 
different polarizations A = 1 and A = 2, as can be seen from (fTHj) . This means that it 
is not necessary to introduce a 2 by 2 matrix analogous to Q as we had to do for the 
quantization of q in (jHJ) (as there is no coupling, this matrix would remain diagonal and 
can therefore be represented by a vector just as well). The Lagrangean of the system is 

|| This is in contrast with vector perturbations, which cannot be generated by scalar-field perturbations 
either. One can easily derive from Q that without a source term the (Oi)-component of the Einstein 
equation gives AS; = 0, so that the vector perturbations are zero. Moreover, even if there is a (non- 
inflation) vector source at some time, one finds from the (zj)-component of the Einstein equation that 
vector perturbations decay as a~ 2 during the expansion of the universe. Hence one only has to consider 
scalar and tensor perturbations in inflation theory, and not vector perturbations. 
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L = \{ip'Ah) 2 ~ |(^ 2 — a " /a){ipAk) 2 and quantization is straightforward: 

kjhiv) = J2~ e i (^Akiv) a\ k + c.c.) , (19) 

A=l 

with the creation and annihilation operators satisfying the usual commutation relations. 
Since the different Fourier modes, as well as the different polarizations, do not couple, 
we drop the subscripts A and k for notational simplicity and consider one generic 
polarization and Fourier mode in the rest of this subsection. 

In a way analogous to (but much simpler than) the treatment for the scalar 
perturbations in subsection 3.2 of ^Hj we can derive the initial conditions for ip and 
the canonical momentum dL/dip' = ip' by using the canonical commutation relations 
between ip and if)' and the condition that the Hamiltonian does not contain any particle 
creation or annihilation terms initially, when k 2 is still much bigger than any other 
scale. This leads to the relations ijj*ip' — ipip'* = i and (ip') 2 + k 2 tf) 2 = 0, with the 
solution ipi = e ia /\/2k and ip' { = i\Jk/2e ia . Here a is an arbitrary phase factor, which 
in a way completely analogous to the scalar case can be shown to be irrelevant to the 
physical correlator, just as the whole sub-horizon region, where ip is simply oscillating, 
is irrelevant. Hence we take a = without loss of generality. 

Realizing that a" /a = 7i 2 (2 — e) we see that the whole treatment in subsection 3.5 
of [TUj of the scalar case is easily applied to this well. In a sufficiently small 

interval around 777^ (the time when k = TC) we find to first order in slow roll with z = krj: 

Hz) = yj^TzB^z) = - m [1 + ~e H (l - 7 - In 2)] (_j , (20) 

where H^ l > is a Hankel function and the expression after the last equals sign is only 
valid for \z\ l.^f On the other hand, the solution in the super-horizon region where 
k ^ Ti, is given by 

r dz' 

i/)(z) =Ca + Da — . (21) 
Jzn a 

Integrating the expression Hi = 7i 2 (l — e), with e = e~u taken constant (first-order slow- 
roll approximation), once with respect to conformal time we find TC(t]) = — [(1 — e^)^] -1 . 
Integrating once more, using TC = (In a)', we find the following first-order approximation 
for a(z) around 77 = 77^: a(z) = an(z / zh)^ 1 ^ 671 ■ Using the identification procedure 
described in jTU] we see that again it is the leading-order term in the expansion in z of 
the Hankel function in the solution for ip in the transition region that turns into the 
dominant solution for ip in the super-horizon region, i.e. the C term in (j21|) . From this 
we derive an expression for C: C = (y2k a-^) _1 [l + e^(l — 7 — In 2)]. Here we omitted 
some unitary factors that are irrelevant to the calculation of the correlator of the tensor 

% Note that the change from the complex time dependence e lkv /V2k in the sub-horizon region, 
satisfying the quantum commutation relation ip*ip' — xpip' = i, to the real time dependence in (121 )[l at 
the end of the transition region, satisfying the classical (H — > 0) relation ip*tp' — tpip'* = 0, is one way 
to see the quantum-to-classical transition taking place. For more information see (351 136j . 
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perturbations. As the D term in (J21|) rapidly decays compared to the C term, we do 
not have to determine D. 

Using we then obtain the final result for hij at later times (i.e. when we can 
neglect the D term) to first order in slow roll for super-horizon modes that crossed the 
horizon during slow-roll inflation: 

/2 2 
^* = T5/5T t 1 + ^ - 7 - ln 2 )] E e 3 *Ak + cc, (22) 

A=l 

where we used the identity k = a^H-^. This result is also valid after inflation, as long 
as the mode k remains super horizon. We see that this expression for h^ is independent 
of time: the super-horizon h^ is simply constant. Of course this could be seen directly 
from the equation of motion for h^ for k 2 negligibly small (i.e. neglecting the Ahy term). 



3. Inflation and the CMB 



In this section (plus the appendix), expressions for the scalar (adiabatic, isocurvature 
and mixing) and tensor amplitudes and spectral indices and their derivatives are derived. 
The results are discussed, paying special attention to multiple-field effects. Moreover, a 
number of important relations between these quantities are derived, which can be used 
as observational tests to determine the presence of multiple fields during inflation. 



3.1. Spectral quantities from inflation 

Having determined the relevant scalar and tensor perturbation quantities at the time of 
recombination (see Q, (fTUj) and (J22J) we can now compute their correlators. We define 
the quantities \5^\ 2 (X denoting adiabatic, isocurvature, mixing or tensor) as 

oh.3 ot.3 ojL3 

| A ad|2 _ ^ /A2 \ I Xiso|2 _ ^ 1^2 \ |A tens | 2 = /% h ij \ 

\°k I — g n 2\ fcad''™' l°fe I — \ ^feiso/W ) \°k I — g^2 \ V h u k /tree ' 



2k 3 f - - \ 

|5- ix | 2 ^ — (($ fciso $ fcad ) ircc + ($ fcad $ fc iso) 4rcc ) • (23) 

Unfortunately there are different conventions in the literature regarding the 
normalization factor. The normalization used here corresponds with |3T1 l3*H], in which 
papers the relation between this amplitude and the observations is explained extensively. 
The normalization factor in, for example, [34j, which was used in the table with 
numerical values in our paper JU] as well, is 9/4 times larger. 

If 1 5£ | 2 depends only weakly on k, the following approximation can be made: 

\^? = K\ 2 {£f\ (24) 

with k$ a certain reference scale. |<5^J 2 and fix are two constants, called the amplitude 
of the CMB power spectrum and the spectral index, respectively. This approximation 
holds good over a wide range of k if hx is close to zero, which is the case for slow-roll 
inflation (see later). Note that originally the spectral index was defined as nx = nx + 1, 
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except for the tensor perturbations where the above definition was indeed used. As this 
was a rather unfortunate source of confusion, newer papers usually adopt the definition 
given here, but to avoid confusion we have added the tilde. 

With these definitions we find for the amplitudes the following expressions, valid 
up to and including first order in slow roll: 



\5 



ad 1 2 
k I 



I rtens 1 2 

\°k I 



1 

6 
400 

IT 



K 2 




50vT 2 




K 2 


H n 


50vr 2 


£n 


K 2 


H n 



50?r 2 l H L 

b(JTC z 



(l-2e H )(l + U£ e U Pe ) 
+ 2B (2l H + f,\ i + 2^e T 2 U Pe + U^ e 5 H U P 
\l-2i n )V e T V e +2BV e T 5 n V t 

- 2~e H )U T Pe V e + 2B (fj^V e + U T Pe 5 H V e ) 
1 + 2{B - l)~e H 



(25) 
(26) 
(27) 
(28) 



Here U Pe , V e and 5 are defined in and (jnj), and B = 2 - 7 - In 2 « 0.7296. Note 
that Ul e U Pe and V^V; can be of order unity and hence the leading-order expressions 
include the (1 + U Pe U Pe ), the VjV e , and the U Pe V e factors, respectively. The reader 
interested in more detail about the derivation is referred to [TTH 12*2"] . 

The spectral indices fix can be calculated from the expressions for \5^\ 2 above: 



n x = 



dhi|5*| 2 



dln|^| 2 dt 



H 



din 1 5 



XI 2 

fc I 



_ _ . (29) 

din A; fc=fco dt H d\nk dt H H n (l - e n ) 

Here we omitted the explicit k = ko from the last two steps, but of course it should be 
applied there as well. In the last step we used di^/d lnfc = (d In fc/dt^) -1 and Tin = k. 
To work out this expression we need the derivatives of Up e and V e with respect to tu'- 



dU T Pe 



dt 



H 



p Upe -2H H fj^ + 



2e 



dV} 



dt H 

Here we used (Q 



H H U Pe 



H 



-H 



Vh 



-5n) - 2H H fi H e 2 , 



(30) 
(31) 



H 1 ) 



H-h(1 — e-yi + §h)Qh- The latter result 
follows from the fact that Q oc rj^ 1 ^ Sn near 777^ (see e.g. equation (63) in jTOj) and using 
the expression for Tt(r]) in the text below (j2*Tj) . The final results to leading (first) order 
in slow roll are 

2l n + ^ + 2ffk<%U Pe + Ul e 5 H U Pe 



FadLo. 



r^iso I l.o. 



I^mixll.o. 



I^tens I l.o. 



1 + UlJJ 



Pe^Pe 



-2 



v e T v e ' 



TJ T V 



-21 



(32) 

(33) 

(34) 
(35) 
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where l.o. means leading order. The result for the adiabatic spectral index was previously 
derived in our paper When comparing (J25Jl - (J28)l and we see that we can 

rewrite the spectral amplitudes to first order in slow roll in terms of the leading-order 
expressions: 



X\2 



|^|f.o.(l + [^e M ]i.o.- J B[n x ]i.o.). (36) 

This means that using we can immediately extend our results for the spectral 
indices to include the next-to-leading-order (i.e. second-order) terms: 



n x = Kli.o. (1 + e H ) + 



dm 



din A; 



-B 



l.o. 



dnx_ 
dink 



(37) 



J l.o. 



The dnx /d Ink are given in the appendix. 



3.2. Discussion 

We continue with a discussion of the results derived in the previous subsection, leading 
to several important conclusions. Some of these were already given in ^U|, but are 
repeated here for completeness' sake. 

We see that the fix contain only slow- roll terms, and thus the \5^\ 2 do indeed 
depend on k very weakly. In other words, the spectrum predicted by slow-roll inflation 
is nearly scale-invariant. Of course if V e = 0, there are no isocurvature perturbations and 
n iso and n mix become meaningless. The derivatives of the spectral indices dn x /d In k 
are second order in slow roll. This might be too small according to recent WMAP data 
[HE], where a dnx/d Ink of the same order of magnitude as fix is claimed, leading to 
a problem for slow- roll models of inflation in general. However, as explained in jS] , this 
result is only obtained when combining WMAP with other data sets, especially Lyman-o; 
data. Hence this result is as yet controversial, with several authors [IUJ EH 113 HH] 
claiming that a constant fix is perfectly consistent with the data. 

The explicit multiple-field terms in the amplitudes and the spectral indices are the 
contributions of the terms Up e and V e , which are absent in the single- field case (setting 
them equal to zero we obtain the well-known single- field results, see e.g. [34J). There are 
no isocurvature and mixing contributions in the single-field case. Since both Up e and V e 
are to a large extent determined by f} 1 - (see their definitions in ()12j) ). we can draw the 
important conclusion that the behaviour of fj 1 - during the last 60 e-folds of inflation is 
crucial in order to determine whether multiple-field effects are important. On the other 
hand, the fact that Up e depends on fj 1 - should not be taken as an indication that it is 
a first-order quantity: as was shown in the example in Up e can be of zeroth-order 
importance (because of the integration interval). 

The fact that the entropy perturbations act as sources for the adiabatic perturbation 
(the Up e terms) naturally leads to correlations between adiabatic and isocurvature 
perturbations (described by the mixing amplitude), as had been realized before (see 
e.g. 53JEHI35])- Note that even if Up e = there is still one other term in the mixing 
amplitude, although merely of first order in slow roll. Only if fj^ vanishes as well, the 
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correlations are completely absent. However, at least in the context of slow roll the 
situation where Up e — while fj^, V e ^ is not possible, because the e\Qja under the 
integral in the definition of Up e cannot change sign. Anyhow, if fj — everywhere 
during the last 60 e-folds, there are certainly no correlations. (The authors of [13] 
studied the two-field case and found the derivative of the angle that parametrizes the 
influence of the second field on the background trajectory to be the relevant parameter. 
In the two-field limit this parameter corresponds with fj- 1 , but our result is valid for an 
arbitrary number of fields. A general discussion of two-field models has recently been 
given in [Hj-) 

The gravitational potential only depends on the total entropy perturbation S, 
independently of the total number of fields and the actual number of independent 
entropy perturbations. In our basis S depends directly on the e 2 component of q only (see 
the expression for V e in f!12j) or (fTTjl ). (Of course the other components of q influence the 
equation of motion for q 2 and cannot be neglected in general.) This is the main reason 
why we only consider this total entropy perturbation and a total isocurvature amplitude 
in this paper. However, for a complete understanding of the CMB power spectrum it 
will probably be necessary to investigate the individual entropy perturbations as well. + 

Using the concept of slow roll for the perturbations the quantity Up e can be 
rewritten in terms of background quantities only, as was discussed in [TU] (see 
equation (69) in that paper). Because slow roll was used in the derivation, that 
expression is in principle not valid at the very end of inflation. (Note that until this 
point slow roll was only used during a small transition region around the time of horizon 
crossing deep within the inflationary era.) If it does indeed give a bad approximation 
for Up e , for example if fj 1 - grows very large, a more careful treatment of the transition 
at the end of inflation is necessary. However, in other cases the contribution to the 
integral near the end of inflation can be negligible, for example if fj 1 - goes sufficiently 
rapidly to zero. In those cases the details at the end of inflation are unimportant for 
the adiabatic amplitude ()25|) . An example of this latter case was discussed in [TU] . 
Unfortunately V e depends very much on the details of the transition at the end of 
inflation, so that for an accurate calculation a model of this transition has to be assumed. 
It also depends on the behaviour of S during radiation and matter domination, as shown 
in the derivation of (JEj), but this aspect is investigated in section HI However, one can 
still draw the conclusion that, if fj 1 - goes to zero at the end of inflation, the isocurvature 
perturbations are expected to be negligible compared with the adiabatic one (neglecting 
possible amplifying mechanisms during the transition and preheating). 

Compared with the scalar amplitudes, an overall factor of 1/e-n is missing in the 

+ For example, one can prove that if there is not only baryonic matter, radiation and cold dark matter, 
but hot dark matter and/or quintessence as well, their individual entropy perturbations enter the 
expression for the Sachs- Wolfe effect (see [H]). An investigation of the individual entropy perturbations 
in a model with photons, baryons, neutrinos, and cold dark matter was given in |49|. The same 
authors investigated the detectibility of these entropy perturbations in upcoming experiments in |5U] 
and concluded that we will have to wait for the Planck satellite |51| for sufficient discriminating power. 
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expression for the tensor amplitude, showing that the tensor contribution to the CMB 
from slow-roll inflation will generically be smaller than the scalar one. Moreover, 
this means that the tensor amplitude depends only on (to leading order), thus 
allowing a direct determination of the inflationary energy scale. Another important 
difference between the tensor and scalar quantities is that the tensor ones do not 
depend on multiple-field terms. As explained in subsection 12.31 this is because the 
tensor perturbations are not generated by the scalar-field perturbations and hence do 
not see the difference between multiple- and single-field inflation. It is this fact that 
allows us to derive, in the following subsection, an important relation that can be used 
as an observational test for the importance of multiple-field effects. 



3.3. Observational tests and consistency relations 

From the expressions for the spectral amplitudes and indices ((|25 |) -(j28 j) and (|32j) - (|3~oT) ) 

we can derive some important relations between these quantities. Note that we only 
use the general expressions for these amplitudes and indices in the derivations, not the 
specific form of, for example, V e . Therefore the results are valid very generally, and do 
not depend on the details of the evolution of the entropy perturbations after inflation 
(as was explained in the derivation of (|17|) . these only change the definition of V e ). 

We define the tensor to scalar ratio r as r = |<5fc ns | 2 /|<5fc d | 2 - Then we find the 
following relation to leading order: 

400 e H 200 n tens rp 200n tens 

r -~i + ul e u Pe -~~i + ul e u Pe * u^-- 1 -—— ( 38 ) 

This expression is a very important result: it gives a relation between observable 
quantities and the length of the vector Up e that encodes the effects of multiple fields. 
Once observations have become good enough to determine the tensor to scalar ratio r 
and the tensor spectral index ntens independently, this relation offers an observational 
test to check whether multiple-field effects are important or not. 

Using (J3T?)) and (}3*T|) we can extend this result to next-to-leading order: 



T _ 200 rn, 

u Pe U Pe — 1 



l-{B-\)n tei3S +Bh ad + {B-l) ' <b/i 



- - , — ,,„> - ,1.1 • \— - ' - 1 1 7 ' \ / 

9 r L ^tens din A; 

This does not change the above conclusion: it is still a relation that allows the multiple 
field quantity Up e U Pe to be determined from observations.* We can also (in principle 
determine the other multiple-field quantities V^V e and Up e V e from the observations: 



v T v = son ntens u T v = ntens , im 

■ ' ! | £tcns 1 2 y' | £ iso 1 2 ' u Pe v e g | ^tens 1 2 j | £mix 1 2 ' ' 



with similar extensions to the next order as in (J39|) . 

* There is one other effect that might iead to corrections to this formula: non-vacuum initial states 
caused by trans-Planckian physics [52| . However, this effect is expected to be small: deviations from 
the vacuum initial state should not be too large, otherwise the particle background dominates over 
the potential energy of the inflaton and there is no (standard) inflation For a recent review of 

trans-Planckian effects see |54j . 
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Next to the above results valid for an arbitrary number of fields, we can find some 
further results in the case of one or two fields only. These are the so-called consistency 
relations, which are based on the fact that in those two cases we have more observational 
quantities than inflationary parameters. In the single- field case the unknowns are H-^, 
e-yt and fj^, while the scalar and tensor amplitudes and indices give four observational 
quantities. Hence we can express one of these in terms of the others: a consistency 
relation. Conventionally this is written as follows (valid to leading order) [213 EH EH] : 

r = - ~7T^ tens (one field), (41) 

y 

a relation that follows immediately from the single-field limit of ([38)1 . 

In the two-field case, there are four more inflationary parameters: fj^, Up e and V e 
(both are vectors, but in the two-field case there is only one non-zero component), and 
S-h (this is a matrix, but in the two-field case there is only one unknown component). 
As there are also four more observational (in principle at least) quantities, there is still 
one consistency relation, which can be written as (again only valid to leading order): 

r = — !r ™ tens i 1 - r mix) ( tw o field s)- (42) 
Here r 2 , is defined as 



|£mix|4 


(u T Pe v e y 




{l + U T Pe U Pe )V?V e 



(43) 



where the second expression is valid to leading order. In the two-field case where 
the vectors Up e and V e have only one non-zero component (denoted by U and V, 
respectively), the V drops out and one has r mix = U 2 /(l + U 2 ), which leads to the 
expression in (|42jl . This two-field consistency relation was also derived in a different 
way in ^21EB]4 Of course the two-field consistency relation reduces to the single-field 
one in the appropriate limit because r m ; x is zero then. In the case of three or more 
fields the number of inflationary parameters increases (with four in the case of three 
fields: one additional component in Up e and V e and two in 5n) without an increase 
in the number of observational quantities (at least when considering only the total 
isocurvature perturbation). Hence there are no consistency relations in those cases. For 
some consistency relations concerning the dn^/d Ink see the appendix. 

To conclude this section: assuming that we will be able to measure the scalar and 
tensor amplitudes and indices with sufficient accuracy in the near future (which of course 
implicitly assumes that they are large enough to be measured), we can then first use 
equation (|3*8j) (or (J2HJ)) to check if multiple-field effects are significant at all. Next we 
can use the consistency relation (|4~2|) to distinguish between the cases of two or more 
fields. Note that an accurate measurement of the isocurvature and mixing quantities is 
only necessary for the second step, not for the first. An accurate measurement of the 

(t In JS] a second consistency relation for the two-field case was derived as well. However, this relation 
is not valid in general, as was also pointed out in , but only if the inflation model satisfies certain 
specific additional conditions (namely that the (22)-component of the matrix Su is equal to 2e^ + wL). 
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tensor amplitude and spectral index is, however, essential for both steps (for a recent 
study on detectability issues see [57]). 



4. Entropy perturbations 

In the derivation of the isocurvature perturbation it was first assumed that the 
total entropy perturbation 5* is constant during radiation and matter domination on 
super- horizon scales, leading to equation (fT2j) for V e . We then considered the possibility 
that S is not constant, finding that (|1()J1 is still correct, but now with a more general 
expression (fTTj) for V e . In this section, we investigate the total entropy perturbation 
during radiation and matter domination to determine its evolution with time. 



4-1- Total entropy perturbation 

We consider a universe filled with an arbitrary number N of energy components, labelled 
by the subscript i. The different components each have a pressure pi and an energy 
density p iy as well as pressure and density perturbations. For the total pressure p = ^2 { Pi 
and the total energy density p = J2i Pi we have 

w = -, c 2 = %, p + 3Hp(l + w) = 0, w = -3H(c 2 -w)(l + w). (44) 
P P 

The first two expressions are definitions (w is called the equation of state parameter and 
c 2 the sound velocity), the third one is the energy- momentum conservation condition 
D^T^q = 0, and the last equation follows by writing w = (p/p — p/p)p/p- For the 
individual components we define analogous quantities: u>, = Pi/pi and c 2 = pij pi (note 
that in general w ^ J2i w i an d c \ 7^ Yl,i c l) an d find 



Pi + 3Hpi(l + Wi) = 3Hd, Wi = -3H(c 2 - Wi 



(1 + ^)-- 
Pi\ 



(45) 



with Ci a measure of the interactions between the different components satisfying 
J2iCi = 0- I n the following we make two assumptions regarding the separate 
components, which for the rest are completely arbitrary: 

(i) All components behave as ideal fluids with a constant Wi. 

(ii) There are no interactions: Cj = for all i. 

From equation (|45|) we see that this automatically means that c 2 = Wi (the square in 
c 2 is just convention; c 2 can be negative). Moreover, a constant Wi also means that 
Spi/Spi = Wi. 

Let us remark briefly on the assumption of no interactions with regard to a real 
model. One can think of the following situation. Of the multiple scalar fields during 
inflation, one has decayed to all the Standard Model particles, while the other fields 
have decayed to various kinds of dark matter. Then there are by construction no 
entropy perturbations between the Standard Model components on super-horizon scales 
and the absence or presence of interactions here is irrelevant. There are only entropy 
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perturbations between the dark matter and Standard Model components and between 
the dark matter components themselves, where the assumption of no interactions seems 
quite plausible. It is usually assumed (see e.g. |HHJ Ell) that including interactions 
will have the effect of wiping out the isocurvature perturbations. However, it will be 
interesting to check this more carefully in the near future. 

To rewrite S (defined as S = \{5p — c 2 5 p) / \p — c 2 s p)) we need some auxiliary results. 
In the first place we have that p = £V WiPi, so that 

s P Y p(i + «0 

where we used (J43J) and ()44|) for p\ and p. Using this result the numerator of S is 
rewritten as follows: 



dp - c 2 s 5p = ^2(w k - c 2 )5p k 

k 

= TT, — r y^pi{l + wi){w k -wi)5p k 



p{l + w)2 



\^2,PkPi{^ + w k )(l + wi)(w k - wi)S M (47) 



k.l 



with 



Ski = 7T~~ 7 , — r. (48) 

p k {l + W k ) pi(l+Wl) 

In the last step of (}4T|) we symmetrized the expression in k and /. Completely analogously 
we find for the denominator 

P-c 2 s P = J2(wk ~ c 2 )p k = — - — ~V»K - wi) 2 . (49) 

Our final result for the total entropy perturbation is then 

a lT,k,iPkPii 1 + w k){ 1 + wi){w k -w l )S k i 
4 2^k,l PkPi{w k - Wif 

The S k i are the individual entropy (or isocurvature) perturbations [SHI EH ESj (see 
also |61J). They are antisymmetric in k and I, and in addition one has S k i = S km — Si m . 
This means that the matrix S k i contains (N — 1) independent elements. One can take 
a single reference component and define S k = S k o, so that Ski — S k — Si, with of 
course = 0. Hence if we have a system of N components, there are in general 1 
adiabatic and (N — 1) entropy perturbations. The combination S of these (N — 1) 
entropy perturbations that enters as the source term into the equation for $ is what we 
call the total entropy perturbation. In the case of inflation with iV fields, the adiabatic 
perturbation corresponds with the e\ direction in our basis, while the perturbations 
in the (N — 1) other directions are isocurvature perturbations. The total entropy 
perturbation S then corresponds exactly with the €2 direction in our basis, see (|15j). The 
entropy perturbations S k i are gauge-invariant by definition in the absence of interactions, 
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see e.g. [60 , and can even be defined in such a way that they are gauge-invariant when 
interactions are included, see |61| . The total entropy perturbation 5* defined in (I14|) is 
always gauge-invariant. 

Working out ()50|) in the case of an arbitrary number of matter components labelled 
by m (e.g. baryons or cold dark matter with w m = 0) and an arbitrary number of 
radiation components labelled by r (e.g. photons or hot dark matter with w r = 1/3), 
we obtain 

q_ Ylm,r PmPrSmr _ PmS m Y2r PrS T . . 

2^m,r PmPr t-im P m P r 

In the last step we have singled out one of the radiation components, for example the 
photons 7, as the reference component, so that S m = S mi and S r = S n . (One could 
just as well choose one of the matter components as reference. The only difference 
with (pTTj) is then an overall minus sign.) In the case of a simple two-component system 
consisting of photons 7 and one cold dark matter component C, which is the case usually 
considered in inflationary literature (see e.g. [SHJ HH 113 E3)> this result simplifies to 
S = S c = {Spc/pc) - f (<WPt)- 

4-2. Time dependence of S 

Next we derive an expression for the time derivative of S. First we need an equation of 
motion for Spi. This equation can be derived by working out the condition D^T^ = 
to first order in the perturbations. However, as we are only interested in the super- 
horizon modes, it is simpler to use the method of varying the background equation. (A 
description of this method can be found in e.g. subsection 3.4 of |fO] and references 
therein.) From ()45j) . with the two assumptions of ideal fluids without interactions, we 
then find 

5 Pi + 3H6 Pi (l + - 3$pi(l + Wi) = 0, (52) 

using 5H = (5 In a)' = — $. Equation (|52*|) can also be found in [HOI- From this result 
together with (|4*3jl we easily derive that 

= fto £SN_ (fc „ = a (53) 



p k (l + w k ) \5p k p ki 

When differentiating S, given in (J5U|) . this means that the time dependence is completely 
determined by the background quantities. We find after a long calculation 

~ 3 ^Yjijki \i Wk + Wl ) ~ ( Wi + w i)K w k - wi)(wi - wj) 2 (l + w k )(l + wi)pipjp k piSki 
S = —H- 
4 



(j2i,j(wi-Wj) 2 PiP3 



~ 4 (l + w)(w-c 2 s ) 2 p 3 ' [ j 

The second form of the result comes about after inserting the definition of S k i (|48|) . 
substantial index manipulation in the numerator and using (|49J) in the denominator. 
It is more compact, but in some ways the first expression is more useful. From the 



Multiple-field inflation and the CMB 



18 



first expression we can see immediately that the time derivative of S will be zero in the 
case of components with only two different values of W{ (one cannot make all three of 
[(wk + wi) — (wi + Wj)], (wi — Wj) and {w k — w{) unequal to zero in that case). Hence 
we can draw the important conclusion that in a universe consisting only of matter and 
radiation components, S remains constant on super-horizon scales, irrespective of how 
many kinds of matter and radiation there are (provided that the two assumptions of 
constant tOj and no interactions are valid). This includes a universe with an arbitrary 
number of hot and cold dark matter components. 

Let us also give the result for S in the case that we relax the second assumption, 
that is, if we include interactions. The sound velocity and total entropy perturbation 
are then given by 

c 2 Pi(l + Wj) -Cj w - _ 1 Efc/ 1 + wi){w k - wi)pi5p k + wiCi5p k 

s p(l + w) " 4 Yjk,l( l + W l)( W k-Wl)PkPl + WiCip k 

(Since S k i is no longer zero, it is not useful to rewrite S in terms of the individual 
entropy perturbations in this case.) After a long calculation with substantial index 
manipulation, during which one should keep in mind that Yli @i = 0, we finally obtain: 

S = -H [(1 +w)(w- c 2 s ) 2 p 3 ] 1 y^|wj(wj - Wj)(wl - WiW k + WiWj)piPj5p k 

+Wj(wi - w k )(l + 2wi - wj + w k )piCjbp k - WiWjCiCjdpk 

+Wj(wi - w k )pi^5p k - WiW k (wi - Wj)piPjbC k + WjWkPiCjdCkj. (56) 

A further treatment of the interacting case is postponed to a future publication. 

Next we consider what happens in the case that there are components with three 
different values of w^, for example matter, radiation and a cosmological constant, again 
neglecting interactions. With more than two different components S is in general not 
constant. Even though (J54)) gives an expression for its time derivative, it is in general 
difficult to solve explicitly because the time dependence of the perturbations 5p k is not 
trivial to determine: from (}52|) we see that it depends on $ so that the equations for 
S and $ have to be solved together. However, it turns out that in the special case of 
components with three different values of Wj, the time dependence of S is determined 
by the background energy densities only and does not depend on $. 

The main idea is that we take another time derivative and calculate S. For three 
components we find that the numerator of S in (jSljl can be written as 

g(wi, w 2 , w 3 )pip 2 5p3 + g{w 3 , w 1 ,w 2 )pip35p 2 + g(w 2 , w 3 , wi)p 2 p 3 Sp 1 , (57) 

with g(a, b, c) = (a — b) 2 (c 2 + ab — c(a-\-b)). Taking the time derivative of this expression 
using f45|) and (J52|) we find simply — 3H(3 + w\ + w 2 + w 3 ) times the same expression 
(|57|) (the $ terms from exactly cancel, which can even be proved to be true for 
an arbitrary number of components). Hence we conclude that in the three-component 
case S oc S, where the proportionality factor only depends on the background energy 
densities, not on the energy density perturbations. This means that, given initial 
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conditions for S and S, one can explicitly solve for S(t), without needing solutions 
for Spi(t) and $(£). For four or more components this is no longer true. 

As an explicit example we consider a system with matter m, radiation r and a 
cosmological constant A, where w m = 0, w r = 1/3, and w\ = —1. Then we have 

q A q P 2 ™Pr + PmPr ~ 9 PmPA ~ ^PrPA ~ 23p m p r p A 
O^N — ^bN —j — — — — r . 

(Pm + lPr){PmPr + 9p m p A + 16p r p A ) 

Even though this is not during inflation, it turns out to be convenient to use the number 
of e-folds A" as time variable to remove the H and its derivative. We choose A" to be zero 
at the present time and negative before that. The functions p-i{N) can be determined 
from Pi(N) = Qip c exp(— 3(1 + Wi)N), where f2« = pj(0)/p c is the present density 
parameter of component i and p c = 3Hq/k 2 is the present critical density, which drops 
out of the equations, however. Using data from WMAP j2j we have 

tt m = 0.3, tt r = 5 • 10" 5 , tt A = 0.7, N eq = -8.7, N rec = -7.0, (59) 

where the subscript 'eq' denotes matter-radiation equality. 

Equation (|58jl can easily be solved numerically. However, using the approximation 
that p r ^> p m ^> pa for N* < N < N eq and p m ^> p r ^> p A for N CCL < A^ < N Tec we can 
also find an analytical solution that agrees very well with the exact numerical one. The 
result is 

S(N) = \ § * + ^ ~ X ) fOT N ^ N *> (m) 

1 ) \S*- 1§* N + JS>e 3 ^--^ (| + e 4 ^-)) for N cq < N < iV rec , 1 ) 

where the superscript * denotes evaluation at the beginning A 7 * ~ —60 of the radiation- 
dominated era. To determine if the time dependence of S leads to significant effects we 
must have an estimate for S* N . Evaluating and (J5U|) at A 7 * (where p r ^> p m ^> p A ) 
for the situation under consideration and combining them we find 

S* N = -48^e 3N *S*. (61) 

Here we have taken d~p\ = 0, assuming A to be a pure cosmological constant. This 
means that 

S(N Tec ) = S* - 12 ^ e ^-N cq S* = (1- 10~ 7 ) S*. (62) 

Hence the effect of the fact that S is not constant in this system is completely negligible. 
This system is practically equivalent to one without a cosmological constant, which was 
of course to be expected as the energy density of the cosmological constant is so much 
smaller than that of matter and radiation before recombination. (As soon as one takes a 
third component with ^>p\j p\ ^> S* or with u> 3 ^ —1, it is no longer possible to express 
S* N in terms of S* only. Even though the solution for S(N) can still be calculated 
explicitly in these cases, it becomes more difficult to determine the relative importance 
of the time dependence of S.) 
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5. Summary and conclusions 

In this paper, we investigated some consequences of multiple-field inflation for the cosmic 
microwave background radiation. Building on the theory of ^| we derived expressions 
for the amplitudes (JUJ-fllU) and the spectral indices (P2M3SI), jSZ|) of the CMB power 
spectrum, all valid to next-to-leading order in slow roll. We also derived expressions valid 
to leading order for the derivatives of the spectral indices (|A.lj) - (jA.4j) in the appendix. 
All this in the context of a very general inflation theory with an arbitrary number of real 
scalar fields that may be the coordinates of a non-trivial field manifold (i.e. have non- 
minimal kinetic terms). There are four different versions of all these spectral quantities: 
three related to the scalar perturbations (adiabatic, (total) isocurvature and mixing 
between those two) and one related to the tensor perturbations. 

These expressions were discussed in subsection lrS.2l To summarize, multiple- 
field effects can be important for the scalar spectral quantities, not only for the 
isocurvature and mixing components (which are absent in the single-field case), but 
also for the adiabatic ones. In all multiple- field terms the slow-roll function rj , which 
measures the size of the acceleration perpendicular to the field velocity, plays a key 
role; if it is negligible during the last 60 e-folds of inflation, multiple-field effects 
are unimportant. Unfortunately, to work out the expressions for the multiple- field 
terms explicitly, especially for the isocurvature perturbations, a careful analysis of the 
transition at the end of inflation as well as of the era of (p)reheating is in general 
required. However, as was shown in [HI], there is a wide class of models where fj ± 
goes to zero at the end of inflation (while being non- negligible before that). Then the 
integral expression for the multiple-field contributions to the adiabatic perturbation can 
be worked out explicitly without knowing the details of the transition, while isocurvature 
perturbations are expected to be unimportant in this case (barring possible amplification 
mechanisms during preheating). Even in those models multiple-field effects in the 
adiabatic component can be of leading-order importance. The leading-order terms of the 
spectral indices are of first order in slow roll, while their derivatives are of second order. 
Hence multiple-field slow-roll inflation generically predicts a CMB power spectrum that 
is close to scale- invariant. The tensor spectral quantities do not depend on multiple- field 
effects. 

From the expressions for the amplitudes and spectral indices we derived some 
important relations, which are valid very generally and do not depend on the details 
of what happens at the end of and after inflation. Most important is equation (|38|) (or 
its extension to next-to- leading order (|39jl). which gives the size of the multiple-field 
contribution to the adiabatic perturbation in terms of the (in principle) observable 
spectral quantities. In other words, this relation provides an observational test to 
determine if multiple-field effects play a significant role during the last 60 e-folds of 
inflation. It does require, however, a sufficiently accurate measurement of the tensor 
spectral quantities. 

In the case of only one or two fields, we have the special situation that there is one 
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more observational quantity than there are inflationary parameters. This means that we 
can then derive a consistency relation between the various spectral quantities, which can 
in principle be checked observationally, allowing one to distinguish between the cases 
of one, two, or more fields, provided that the observations are sufficiently accurate. We 
derived that consistency relation in subsection ECU If the derivative of the spectral index 
can also be measured, there is an additional consistency relation, which was derived in 
the appendix. With all these relations one should keep in mind that, while a large effect 
is probably a clear proof of multiple- field effects being important, a small effect might 
also signify the presence of trans-Planckian or short-distance physics (see e.g. jSSEZI)- 

An important ingredient of the calculation of the isocurvature amplitude was the 
observation that, although there may be many individual isocurvature perturbations 
in a multi-component system, only the total entropy perturbation S enters into the 
equation of motion for the gravitational potential. With the assumptions that the 
various components behave as ideal fluids and have no interactions, the time derivative 
of S on super-horizon scales during radiation and matter domination was worked out in 
([54)1 . It was found that, if there are only two different types of energy in the universe 
(i.e. two different equations of state, e.g. baryons and cold dark matter with p = and 
photons and hot dark matter with p = |p), S is simply constant. Although this is 
no longer true if there is a third type of energy, for example a cosmological constant 
with p = — p, we showed that for the general three-component case the equation of 
motion for S(t) only depends on the background energy densities and can be computed 
explicitly, independently of the gravitational potential. An explicit calculation of a three- 
component case with matter, radiation, and a pure cosmological constant (with realistic 
values of the parameters) showed that the time evolution of 5* before recombination can 
be completely neglected in that case. 

To be more general we also derived, at the end of subsection 12.21 an expression for 
the isocurvature perturbation in the case that S is not constant. It turned out that the 
effects of this can be absorbed in the definition of the multiple-field term, so that the 
general form of the amplitudes and spectral indices is unchanged. This means that the 
observational relations in subsection 13 .31 do not depend on whether S is constant or not. 
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Appendix A. Expressions for dn/d lnfc 

The recent WMAP results have made it interesting to look at the derivatives of the 
spectral indices with respect to k as well. In a straightforward calculation (using (J7J) for 
the derivatives of the slow-roll functions), analogous to the derivation of the spectral 
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indices themselves in subsection 13. 1[ we find to leading (i.e. second) order in slow roll: 
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In these equations the derivative 5n = d5n/dtn occurs. In specific models it can be 
calculated from the potential and the field metric using the definitions ([TTf : 
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The derivatives of the basis vectors are given in equation (|4*|). 

If we also include the dhx /d In k as observational quantities, there are, in the single- 
field case, two of those and one additional inflationary parameter (|fj, so that there 
should be one more consistency relation (in addition to ()41)l ). It can be written as 
1 dn ten 
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In the two-field case, there are two more observational quantities but also two more 
inflationary parameters (£ 2 « and 5n), so that we also have a single additional consistency 
relation: 
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where r 2 iY is defined in ()43|) . 
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